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Low temperature NMR in La 2 Ni0 4 i7 
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The experimental observations of stripes in superconducting cuprates and insulating nickelates 
clearly show the modulation in charge and spin density. However, these have proven to be rather 
insensitive to the harmonic structure and (site or bond) ordering. Using 139 La NMR in L^NiC^+a 
with S = 0.17, we show that in the 1/3 hole doped nickelate below the freezing temperature the 
stripes are strongly solitonic and site ordered with Ni 3+ ions carrying 5 — 1/2 in the domain walls 
and Ni 2+ ions with 5 = 1 in the domains. 

PACS numbers: 76.60.-k, 74.72.Dn, 75.30.Ds, 75.40. Gb 



Stripe phases have by now been observed in a variety 
of doped Mott-insulators, like cuprates, nickelates and 
manganites Q. Nevertheless there is still a remarkable 
lack of knowledge on the details of the charge- and spin 
distributions in these novel electronic phases. At the 
same time there is a growing body of theoretical liter- 
ature dealing with the microscopic origin of stripe for- 
mation, predicting stripes starting from rather different 
physical perspectives [0-(5). Given its potential connec- 
tions for instance to the mechanism of superconductivity, 
it is a matter of high urgency to find out how the stripes 
look like in detail. In this regard the various theories 
lead to quite different predictions. Because the spin- and 
charge distributions are inhomogeneous, NMR with its 
microscopic sensitivity could in principle yield important 
information. However, in cuprates attempts in this direc- 
tion have been frustrated due to the anomalous, glassy 
ordering dynamics obscuring the static signal. Neutron 
scattering shows that both in cuprates and nickelates the 
correlation length of the stripe order is relatively small. 
This disorderly nature comes to play in NMR |||7]] in 
the form of a peculiar loss of signal intensity (wipe-out) 
|]6|-|lT| explained by a spread in very short spin-dephasing 
times T2 [|6| Jlc| , pT|l . Even at the lowest temperatures the 
signal recovery is only partial indicating that this dy- 
namics is still at work at temperatures as low as 0.3 K 
in the cuprates, prohibiting attempts to deduce infor- 
mation about the static order from the NMR data 0. 
We demonstrate here that in a nickelate stripe system 
the nature of the stripe order (established by neutron 
scattering) can be deduced in detail from NMR, despite 
the strong similarities with the cuprate NMR at higher 
temperatures. We find the stripe structure to be strongly 
solitonic, with sharply defined charge stripes with a width 
which is not exceeding the lattice constant by much. Sur- 
prisingly, they look quite like the site centered stripes 



predicted by early mean-field calculations for the nicke- 
late system Q . 

Before analyzing the low temperature line shape, which 
is the central issue of this Letter, we first introduce the 
main features of the NMR line shape, derive the value of 
the hyperfine coupling and explain the partial recovery 
of the signal intensity from the relaxation data. 

The NMR measurements were performed on the same 
La2NiC>4.i7 crystals as measured before |]|l|], by sweep- 
ing the external field B at various frequencies. The crys- 
tal symmetry is almost tetragonal JlJ] , which makes the 
139 La-spectra (/ = 7/2) strongly dependent on the direc- 
tion of the magnetic field with respect to the crystallo- 
graphic axes. 
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FIG. 1. Ordering of excess oxygen, holes and spins. Open 
circles refer to the ordered excess oxygen sites at z = 0.217c, 
which introduce holes assumed to be Ni 3+ with 5=1/2 
(gray circles, site-order scenario). The arrows refer to the 
Ni 2+ spins with 5=1. Black circles indicate the 139 La sites 
at z = 0.14c. The dotted lines are the in-plane boundaries 
for the unit cell before doping. The drawn tilted rectangle is 
the enlarged unit cell after doping. 
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In La2NiO/t.i8, with hole doping very close to 
La2NiO4.1v, details about the excess oxygen atoms po- 
sitions are available [[uj. The interstitial oxygen sites 
are at (0.183a, 0.183b, 0.217c) or equivalent positions. A 
scenario for the excess oxygen ordering is shown in Fig|l]. 
There is one excess oxygen atom for each six unit cells. 
The number of 139 La sites, which are far from the excess 
oxygen atoms, is twice as large as the number closest to 
the excess oxygen atoms. For La 2 M04.i7, using simple 
point charge calculations, we expect the interstitial oxy- 
gen to change the electrical field gradient (EFG) around 
the La-sites and to split the 139 La NMR line into two: 
an A-line due to A-sites with the main component of 
the EFG (V zz ) along the c-axis and a B-line due to two 
times less abundant B-sites with the crystal field gradient 
in the (ab)-plane. This prediction is in agreement with 
the experimental intensity ratio of line A to B Klffq] . 
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FIG. 2. (a) 1 for line A (see text) for m=3/2 at 
44.196 MHz, obtained from stretched exponential fits to the 
recoveries. The dotted line is a fit with Eo = 90 K. Drawn 
line is a fit using a half Gaussian distribution in activation 
energy, (b) T^ 1 versus 1/T at low temperatures. An acti- 
vation energy of 2.3 K fits the data well. Sites with lower 
activation energies will not be observable (m =1/2 refers to 
the +1/2, -1/2 and m = 3/2 to the 3/2, 1/2 transition. 

Above 50 K, the 139 La relaxation rates can be fitted 
by an Arrhenius expression T^ 1 cx exp(Eo/T), see Fig. 
||a [|| . Fits allow for a (full or half gaussian) distribution 
in activation energies around Eq = 90 K with a width 
of A = 25 K. Above the magnetic freezing temperature 
the rates for line A are at least twice the rates of line B, 
and have almost the same temperature dependence ||. 
To extract the amplitude of the fluctuating fields and its 
correlation time from the relaxation data, the value of 
the hyperfine constant is needed. This constant can be 
determined in several ways. In doped La2CuC>4 the static 
internal magnetic field amounts to about 0.1 T |l, which 
is similar to the value found in the undoped case [ p^[ . 
Because the hyperfine field at the lanthanum site in the 
nickelates is calculated to be about 18 times larger than 
in the cuprates p7y , its value will be around 1.8 T. Yoshi- 
nari et al. |1£| reported a value of about 2.3 ± 0.5 T//ig 
for the hyperfine coupling in Sr-doped nickelates with an 



ordered moment of ~ 0.7 ilb- The so obtained internal 
field of 1.6 T agrees well with the internal field of about 
1.78 T reported by Wada et al. @ for La 2 Ni0 4 .io. Also 
by comparing our susceptibility and linewidth data at 
high temperature we can estimate the hyperfine field to 
be about 2.0 T, consistent with the earlier quoted values 
p3[ . Using a hyperfine field of 1.8 T, the fit to the re- 
laxation data gives a value of 2600 s _1 for 7 2 /i 2 >t oc , and 
hence = 1.8 x 10" 11 s. T^ 1 has almost the same T 
dependence as . It deviates only for T > 200 K. 

The relaxation data below 20 K can not be fitted by 
the same parameters, see Fig. ||b. The rates can be 
reproduced by T^ 1 oc 7 2 /i 2 ,exp(-£t)/T)/(u; 2 T 00 ), with 
Eq = 2.3 ± 0.2 K and ho the strength of the fluctuating 
field. Measurements of the rates were performed for dif- 
ferent satellites at 38.2 MHz. Using the same methods as 
described in Ref . [|| , a magnetic origin for the fluctuating 
fields is found at 160 K, 90 K and 1.6 K. 

The wipe-out of the 139 La NMR signal intensity and its 
partial reappearance below ~ 11 K || can be explained 
in the same terms as used before in the cuprates Q in- 
dicating a distribution in spin dephasing times, hence 
activation energies fTo] ] . At low temperatures only nuclei 
relaxing with activation energies > Eq lead to the mea- 
sured signal, while the others will have dephasing rates 
higher than 10 4 Hz down to very low temperatures. Even 
at T sa 0.5 K slow spin fluctuations apparently prevent 
the full recovery of the signal. The intensity recoveries 
can be fitted with A = 2.8 K and E ^ K These 
results are consistent with the relaxation data because 
only nuclei with the higher activation energies will be 
visible. 

What are the implications of these observations? In 
the renormalized classical limit below charge ordering the 
spin correlation rate r(T) will be proportional to the spin 
correlation length £(T) divided by the spin wave veloc- 
ity c: Tf 1 - (£/c)(T/27rp) 2 (l + T/2np] T 2 , with p the 
spin stiffness Ijjj- In La2NiC>4.i7 according to diffrac- 
tion studies the stripes have a finite and temperature 
independent correlation length [po|]2l| l. If this has a dy- 
namic origin, it will be consistent with a distribution in 
activation energies or r's. Possible sources for such a dis- 
tribution in the spin dynamics are dislocations ]22] , ^3]| in 
the domain pattern or elastic deformations J23| . 

We now come to the central part of this paper. We will 
show that although the interstitial oxygens have a direct 
influence on the La signal via the electric field gradients, 
it is still possible using 139 La NMR to probe the charge 
and spin ordering in the Ni02 plane. 

Below 15 K, when the signal starts to regain some of its 
intensity, the NMR lines are largely broadened. In Figs. 

the spectra are shown for B || c and B ± c . Simula- 
tions of the spectra taking into account the quadrupolar 
and Zeeman contributions p4| ] are shown in the lower 
panels. The spectra can be decomposed into 2 La sites, 
one with resolved satellites and the other with overlap- 
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ping lines due to much stronger magnetic broadening. 
Magnetic freezing or ordering is also visible in the bulk 
magnetic susceptibility, which starts to deviate from the 
Curie law and peaks at T = 17 ±3 K The peak 

in —dx/dT indicates a series of cusps at different spin 
freezing temperatures close to 17 K and can be due to 
spin frustration or spin clusters pa] . 
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FIG. 3. High (a) and low (b) T spectra for B \\ c and their 
simulations (lower panels). The simulations give precise val- 
ues for the EFG parameters but are not very sensitive for 
the strength of the internal field (at low T). Parameters are 
discussed in the text. 

We first show that contrary to the observation in the 
cuprates , the low temperature line pattern in the nick- 
elate is not due to motional narrowing. In case of two 
lines separated by Su the condition for motional narrow- 
ing is that the correlation time r for the (spin) 
fluctuations is much less than l/Sto. The splitting Suj 
is given by 7„ times the value of the internal field ho'. 
Suj = 7„ < ho >. The value of r is coupled to the relax- 
ation rate via r ~ 7 2 < h\ > Ti/u 2 . Hence the condition 
can be rewritten as R m = 7^ < h\ > 3 / 2 Ti/lj 2 < 1. Us- 
ing typical values for ho = 1.5 T, u> — 2tt x 40 MHz and 
Ti = 50 us, R m ~ 10 2 . 

What is the information we might extract from the 
lineshape data? The high T data will be sensitive for the 
EFG parameters and the low T data in addition for the 
internal or local field. Since the internal field lies in the 
a&-plane, the spectra for B \\ c will be rather insensitive 
to the local field. For B \\ ab the antiferromagnetic align- 
ment of the electron spins we expect to lead to a splitting 
of the lines. The simulations of the spectra are shown in 
the lower panels. The principal axes of the EFG with re- 
spect to the crystal axes are described by Eulerian angles 
(a, /3, 7) |^6| and the external magnetic field is described 
by polar angles {9, <p) with respect to crystallographic 
axes. 
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FIG. 4. High (a) and low (b) T spectra for B \\ ab and 
their simulations (lower panels). Spectra are simulated with 
the same parameters as in Fig.H. The internal field leads to a 
splitting of the lines and hence can be determined precisely. 

For B || c (0 = 0) the high T data (e.g. at 75 K) 
at v = 38.2 MHz are well reproduced by a quadrupolar 
splitting vq = eQV zz /4hI(2I — 1) (with Q the electric 
quadrupole moment of the 139 La-nucleus) of 4.5 MHz. 
The anisotropy parameter rj = (V xx — V yy )/V zz = 0. 
The width of the satellites in the spectra can be sim- 
ulated by introducing a spread in the EFG parameter 
A^q = 0.5 MHz, while the width of the main line requires 
a (dipolar) field distribution of 0.05 T at the La(A)- 
site. For the La(B)-site these values are vq=8.5 MHz, 
t] = 0.75, Aisq = 1.0 MHz and A B = 0.05 T. The other 
parameters are for the A-site: a = 0, (3 = 0, 7 = 0; 
for the B-site: a = 0, j3 = ir/2, 7 = — ir/6. According 
to the simulation of the low temperature data the main 
effect of the magnetic freezing is the increase of the mag- 
netic broadening (As) to 0.5, resp. 0.15 T. The spectra 
shown use internal fields of 1.5 T at the A-site and 0.5 T 
at the B-site. These internal field values follow from the 
simulations for B _L c, where the values are much more 
restricted. For the B _L c orientation simulation of the 
f07 K data with the same parameter set as for B \\ c 
reproduce the data well (ip = t/6). For the 4.2 K data 
the antiparallel alignment of the internal field leads to 
a splitting of the lines and is described by polar angles 
(6 = ir/2,(p = 7r/3±7r/2) with respect to crystallographic 
axes. 

To explain these remarkable differences, let us assume 
that we deal with site centered stripes. These stripes 
naturally lead to two kinds of Ni-ions: Ni 2+ (A-sites) 
and Ni 3+ (B-sites) with two kinds of spins 5=1 and 
5=1/2 respectively. As a consequence the La sites will 
experience different hyperfine fields. The internal dipolar 
magnetic field arising from the Ni-spins at the La-sites 
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will be of the order of 0.2 T. However, the La nuclei just 
above and below Ni 2+ in addition will have an exchange 
coupling via the oxygens. This hyperfine field due to the 
overlap of the Ni 3d 2 2 and La 6s orbitals through the 2p z 
orbital of apical oxygen, is about 1.8 T in the undoped 
samples and is almost doping independent \27\. Due to 
the different occupation of the 2>d z 2 orbital of Ni 3+ the 
exchange coupling between the 5 = 1/2 Ni 3+ spins and 
the La sites will be weaker J28j,^9|. The ratio between 
139 La sites that feel the exchange field of Ni 3+ and Ni 2+ 
will be close to 1:2. This difference in intensity ratio 
and hyperfine field are indeed the main characteristics of 
the line shape and hence are well accounted for by this 
scenario. Can bond centered stripes explain the obser- 
vations as well? The hyperfine fields will have the same 
maximum, but the distribution will be different. The line 
shapes for B _L c puts a limit to the field on the B sites 
of at most 0.3 T, which rules out this possibility. 

The experiments show that apart from the internal 
field we need to introduce an appreciable magnetic broad- 
ening. Part of the broadening might be due to canting 
of the spins in the ordered phase in the Ni02-plane away 



from the charge and spin stripe direction |21|, which ef- 
fect we have not included into this calculation. Another 
reason for the extra broadening might be found in the 
finite size of the correlated magnetic regions. 

Summarizing, in La2Ni04.i7 interstitial oxygens deter- 
mine the line profile above the wipe-out temperature. 
The NMR intensity loss above the spin freezing or or- 
dering temperature around 20 K is linked to a spread in 
spin-dephasing as in the cuprates. From the angular and 
temperature dependence of the La line profiles, we show 
that the distribution of the internal fields is in agree- 
ment with two kinds of Ni-sites with different ionicity 
and hence different hyperfine interaction with the visible 
La sites. Site centered stripes of the kind predicted by 
mean- field theory jl2] fit the low temperature data of the 
visible La-nuclei remarkably well. 
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